function [x,k]=Gaussmethod(~,~,~,~,~)
%A是线性方程组的左端矩阵
% b 是右端向量
%x0 是迭代初始值
% N表选代次数上限，若选代次数大于N，则选代失败% % emg 表示控制精度
用Gauss-Seidel迭代法求线性方程组Ax=b的解%k表选代次数
%x表示用迭代法求得的线性方程组的近似解
n=length(A);
xl=zeros(n，1)；x2=zeros(n,1)；
x1=x0；
r=max(abs(b-A×xl))；
k=0；
while r>emg
for i=l:n
sum=0；
for j=1:n
if j>i
sum=sutn+A(i，j)*xl(j)；
elseif j<i
sum=sum+A(ij)*x2(j)；
end
end
x2(i)=(b(i)-sum)/A(i，i)；
end
r=max(abs(x2-xl))；
xl=x2；
k=k+1；
if k>N
disp(迭代失败，返回)；
return;
end
end
x=x1; ，_4.+r=1